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A semiempirical shell model mass equation based on magic number Z = 126 and 
applicable to translead elements with TV < 126 is presented. For a-decay energies the 
J> ■ equation is shown to have a high predictive power and an rms deviation from the data 
of about 100 keV. The rms deviations for masses and other mass differences are between 
O ■ about 200 and 300 keV. 
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Recent progress in superheavy elements (SHE) research reaching to 293 118 and its a- 
decay products |J makes it necessary to find an appropriate substitute for the semiempiri- 
cal shell-model mass equation (SSME) || (see also ref. |§) for nuclei in the neighbourhood 
of Z = 114 and beyond |4j]. The a energies of the decaying chain vary smoothly from 
293 118 to 269 Sg (Z = 106), with no indication of magicity at Z = 114 in these nuclei [|TJ, 
whereas the SSME assumes that Z = 114 is the next spherical proton magic number after 
lead and it stops there. Furthermore, the SSME becomes unsuitable for extrapolation 
already earlier, beyond Hs (Z = 108), as shown by its increasing deviations from the data 
when Z increases. (Like in fig. 4 of ref. 

Recent phenomenological studies of BE(2) systematics f| and of the persistence of the 
Wigner term in masses of heavy nuclei indicate Z = 126 as the next spherical proton 
magic number after lead, and this is consistent with considerations based on nuclear dif- 
fuseness [|7| . Recent self-consistent and relativistic mean-field calculations PJT2| variously 
predict proton magicities for Z = 114, 120, 124 and 126. 



During the early stages of the SSME [jTSJ, when it was adjusted separately in individual 
shell regions in the N — Z plane, both Z = 114 and Z = 126, which were at the time 
considered possible candidates for the postlead proton magic number (see, e.g., ref. jl4]|), 
were tried as a shell region boundary in each of the two heaviest regions with Z > 82 and 
respective iV-boundaries 82 < iV < 126 (called here region A) and 126 < N < 184 (called 
region B). The agreement with the data was about the same for both choices, and the 
prevailing view in the mid nineteen-seventies led to the choice of Z = 114 for the SSME 
mass table [[J. In a recent communication ||J] we showed that the early Z = 126 results 
have a high predictive power in the interior of region B and proposed their use there as 
a predictive tool in SHE research. In the present note we study the predictive power or 
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extrapolatability of the early Z = 126 results in region A and propose a mass equation 
based on them as a substitute for the SSME [Q] in the interior of the region. The study is 
based on the newer data that became available after the adjustments were made, as was 
done in refs. I5[-|TB| . 



In the SSME the total nuclear energy E is a sum of pairing, deformation and Coulomb 
energies: 

E (N, Z) = E pavr (N, Z) + E def (N, Z) + E Coul (N, Z) . (1) 
The form of Ec ou i is the same in all shell regions: 

f17 \ ^3 

E Coul (N, Z) = \-J>) [a c + P C (Z- Z ) +! C {Z- Z ) 2 } , (2) 

and that of E pa i T is the same separately in all diagonal shell regions, where the major 
valence shells are the same for neutrons and for protons, and in all non-diagonal regions, 
where the neutron and proton valence shells are different. Unlike in 0, with Z = 126 
rather than 114 as an upper proton boundary, region A becomes a diagonal region with 
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(N, Z) = 
a + (3 {A - Aq) + 7 {A 
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The part E de f for region A with Z = 126 as upper proton boundary is 
E def (N, Z) = {^j [<p n $ n (N, Z) + ^ 2 o [^20 (N, Z) + ^ 20 (Z, N)]] 



with 



$ n (N, Z) = {N- 82) (126 - N) {Z - 82) (126 - Z) 
^ 2 o (N, Z) = (N - 82) 2 (126 - N) 2 (N - 104) . 



(4) 

(5) 
(6) 



In eqs. (@)— fl|) A = N + Z and T = \T Z \ = ~\N — Z\f\. The respective values of Z and 
Aq are 82 and 164. The coefficients multiplying the functions of iV and Z are adjustable 
parameters which were determined by a least-squares adjustment to the data Their 
numerical values resulting from that adjustment are given in the second column of table 
I. The mass excesses AM(N, Z) are obtained by adding to eq. (|Tj) the sum of nucleon 
mass excesses NAM n + ZAM#. 

The experimental data used in the adjustment included 29 masses and 62 Q a values 
connecting unknown masses (Ref. p0| augmented by data from the literature up to Spring 
1973). Presently there are 150 known experimental masses in region A and 3 Q a values 



2 In the as yet unknown odd-odd N = Z translead nuclei the g.s. is expected to have T = \T Z \ 
seniority zero, whereas eq. (|l|) with T = \T Z \ gives the energy of a low excited seniority two state 
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connecting unknown masses (Refs. (excluding values denoted "systematics" (#)) and 
[p2| , augmented by data from the recent literature). There are 121 new masses that were 
not used in the adjustments. 

Fig. 1 shows the deviations from the data of the predictions of eq. ([[]) (with the 
definitions @ — (§) and the coefficients from the second column of table I) for all the 150 
known masses. The deviations are plotted as function of the distance from the line of /in- 
stability, denoted "neutrons from stability" (NFS) and defined by NFS = N — Z — ^+200 
llB"] . Empty circles denote the deviations of the 29 originally adjusted masses and full 
circles mark the deviations of the 121 new masses. The original deviations are relatively 
small. The new deviations are mostly considerably larger and are almost all negative. 
There is a large scatter of the points, superposed on an overall decreasing trend when NFS 
increases in the negative direction. For closer scrutiny we show in fig. 2 lines connecting 
deviations of Pb, Bi, Po and At isotopes as function of T rather than NFS. Isotopic lines 
of other elements follow closely. There is an odd-even staggering of the points, where 
for even-even and odd-odd nuclei the deviations are respectively higher and lower than 
for odd-mass ones, and there is an overall increasing overbinding when T decreases away 
from stability. The increasing staggering when T decreases as compared to the old data 
indicates increasing pairing parameters (see also ref. ||23|| ) and (even more so) k towards 
the proton drip line. Similarly, the increasing binding when T decreases indicates that 
the symmetry-energy parameter e decreases towards the drip line. 

Table II, patterned after similar more elaborate ones |16| , |17| (see also ref. J3|), shows 
the values of 5 av and 5 rmsi the respective average and rms deviations of eq. ([[]) from the 
data, for AM, S n , S p , Qp- and Q a . The deviations are shown separately for the older data 
that were used in the adjustment and for the newer data. The last column shows the error 
ratios S rms : S rms . 

For the older data the 5 av are few tens keV at most, and the 5 rms are in the range 
100-250 keV. For the newer nuclei the agreement with the data of predicted AM (fig. 
1), S n , S p and Qp- values has much deteriorated. On the other hand, the S rms of the 
predicted Q a values has even very slightly improved. The Q a deviations, both old and 
new, are remarkably small. 

The high degree of extrapolatability of the Q a values as compared to the poorer pre- 
dictions of AM, S n , S p and Qp~ is presumably due to the composition of the old data 
set used in the adjustments. The coefficients (3, 7, /3 C , 7 C , (p u and 1P20, which contribute to 
both mass and Q a values, and also a c , which was determined using (3 C and 7 C |L3| (see also 
ref. 0]), were determined by all the 91 available mass and Q a data. On the other hand, 
neglecting their ^-dependence the coefficients a, e, O and k cancel in Q a and they were 
determined essentially by the smaller group of 29 masses found in the nearest-to-stability 
corner of the shell region, with too-small values of O and k and too large value of e as 
compared to more proton- rich nuclei. These values, which to a large extent cancel in Q a , 
are responsible for the large deviations of the predicted new masses. 

One would like to restore to the new AM, S n , S p and Qp- predictions the same degree 
of agreement with the data as the old predictions had, while at the same time retaining 
the high quality of Q a predictions. For eq. (P this goal might most simply be approached 
by making a least-squares adjustment of eq. (]!]) to all the 150 known masses, with only 
four adjustable parameters a, e, O and k, while the other seven coefficients are held fixed 
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on their old values from the second column of table I. The re-adjusted parameters a and 
e would correct the systematic overbinding when T decreases, and the re-adjusted and 
k would decrease the odd-even staggering of the deviations. 

The re-adjusted values of the coefficients a, e, O and k are given in the third column of 
table I. Their changes compared to the second column are in the anticipated directions. 
The respective average and rms mass deviations obtained in the adjustment for all the 
150 mass data are 2 and 246 keV. 

Fig. 3 shows the deviations from the data of the predicted mass values resulting from 
the new adjustment, similarly to fig. 1. For ease of comparison the same nuclei are 
denoted by the same kind of circle (empty or full) in the two figures. The deviations for 
the post 1973 measured nuclei (denoted by full circles) and their odd-even staggering are 
considerably smaller than in fig. 1, and the overall decreasing trend when T decreases has 
largely disappeared, as was expected. The deviations for the older data have worsened, 
though, and there is a new overall oscillatory trend not observed in fig. 1, indicating that 
enlarging eq. ([I]) by additional particle- hole symmetric Edef terms (eq. (PS)) might reduce 
the deviations further. 

Table III shows the resulting 5 av and 5 rms values of the predicted deviations. Like in 
table II they are shown separately for the old 1973 data nearer to stability and for the 
new data extending into the interior of region A towards the proton drip line. For S n , S p 
and Q a the new deviations are similar to those of the original nuclei in table II. For AM 
and Qp- the 5 rms are respectively 1.5 and 1.1 times larger. 

The new mass deviations shown in table III are about one half to two thirds of the 



corresponding deviations of several recent mass models [16|-[T8],^|26[ , and for Q a they are 
smaller. The smaller values of the deviations in table III are presumably due mainly to 
the inclusion of the particle- hole symmetric configuration-interaction terms Ed e f (eq. (fD) 
in eq. (|I|) H. Until a new adjustment of the Z = 126 SSME to the data in both regions 
A and B is undertaken, which would further reduce the deviations in both regions, we 
propose the use of eq. (H), with the new values of a, £, and k and the old values of the 
other coefficients, given in table I, as a predictive tool for masses and their differences , 
and particularly for the extrapolatable-proven Q a values, in the interior of region A. 

It should be emphasized, though, that the above results are not a proof of superior 
magicity of Z = 126 in region A as compared to other recently proposed prediction 
[PHHT2[], because no comparative studies of this kind were made. (See also ref. ||.) 

We thank the Atomic Mass Data Center in Orsay for ref. |22| and Stelian Gelberg and 
Dietmar Kolb for help with the calculations. 
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Table 1 

Values of the coefficients of eq. (f ) as determined by adjustment to the old data [13] (Old 
Value) and as readjusted in the present work (New Value). 



coefficient Old Value (keV) New Value (keV) 

-1.987628 x 10 6 



4.585496 x 10 2 
1.2183 x 10 3 
2.1937 x 10 3 



a 


-1.9902575 x 10 6 




-2.4773664 x 10 4 


1 


-8.51085 x 10 1 


e 


4.658516 x 10 2 





9.762 x 10 2 


K 


1.4965 x 10 3 


of 


7.968418 x 10 5 




2.032906 x 10 4 


7 C 


9.819137 x 10 1 




-4.794 x 10~ 2 


^20 


9.095 x 10~ 4 



Table 2 

Numbers of data N, average deviations 8 av , and rms deviations 5 rms , for eq. (1) with 
the old values of the coefficients from table I. The last column shows the error ratios 

Knew . sold 
"rms ' rms' 



Original nuclei (1973) New nuclei (1973-2000) 







fiav 


^rms 




fiav 


firms 


Error 


Data 


N 


(keV) 


(keV) 


N 


(keV) 


(keV) 


ratio 


AM 


29 


-29 


146 


121 


-807 


1008 


6.88 


S n 


18 


39 


214 


120 


-120 


406 


1.90 


Sp 


22 


9 


182 


104 


132 


417 


2.29 


Qf3- 


15 


-20 


242 


101 


248 


583 


2.41 


Qa 


78 


5 


103 


31 


40 


89 


0.87 



Table 3 

Numbers of data N, average deviations S av , and rms deviations fi rm s, for eq. (1) with the 
new values of the coefficients a,e,Q, k and the old values of the other seven coefficients 
from Table I. 



Original nuclei (1973) New nuclei (1973-2000) 







fiav 


firms 




fiav 


firms 


Data 


N 


(keV) 


(keV) 


N 


(keV) 


(keV) 


AM 


29 


-193 


344 


121 


48 


216 


S n 


18 


158 


416 


120 


-10 


205 


S p 


22 


-144 


202 


104 


18 


184 


Qp- 


15 


-257 


475 


101 


15 


277 


Qa 


78 


-5 


104 


31 


18 


85 
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Figure 1. Deviations of the predicted masses [13] from the presently known data in region 
A. Shown as function of Neutrons From Stability (NFS). 
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Figure 2. Isotopic lines of predicted mass deviations [13] of Pb, Bi, Po and At nuclei 
measured after the adjustments were made. Shown as function of T. 



8 



> 
0) 



Q. 

X 
0) 

< 
I 

o 

(0 

o 



\ i r 







■1 - 



J L 



1 I I I I | I I I I | I I I I | I n 



° % 

Ax 

<5> 



<D 



O 



o 
o 
o 



o 



o 
o 



o 



o 



o 



1 



J I I L 



o O 

o mass known in 1973 o o> 
• mass measured after 1 973 

i i i i i i i i i i i_ 



o — 
J L 



15 -10 -5 

Neutrons from (3-stability NFS 



Figure 3. Deviations of the predicted masses resulting from the present adjustment from 
the presently known data in region A. Shown as function of NFS. 



